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Abstract 

The main constraint on relative position of geometric objects, used in spatial planning 
for computing the C-space maps (for example, in robotics, CAD, and packaging), is the 
relative non-overlapping of objects. This is the simplest constraint in which the minimum 
translational distance between objects is greater than zero, or more generally, than some 
positive value. We present a technique, based on the Minkowski operations, for generating 
the translational C-space maps for spatial planning with more general and more complex 
constraints on the relative position of geometric objects, such as constraints on various types 
(not only on the minimum) of the translational distances between objects. The developed 
technique can also be used, respectively, for spatial planning with constraints on transla- 
tional distances in a given direction, and rotational distances between geometric objects, as 
well as for spatial planning with given dynamic geometric situation of moving objects. 

Keywords: Spatial planning, Configuration space, Minkowski operations. 
1 Introduction 

Problems concerning the relative placement of geometric objects are called spatial planning 
problems [29]. Such problems are important in robotics [24], collision detection [21], [27], and 
computer-aided design and manufacturing (CAD/CAM) [11], [54]. 

A technique commonly used for solving spatial planning problems is the configuration 
space (or C-space) approach, based on representing each placement of an object, i.e., its posi- 
tion and orientation, as a point in some parametric C-space [28], [29]. (Each coordinate of the 
C-space represents a degree of freedom in the position or orientation of the object.) 

Given a collection of objects, the translational spatial planning problem consists in comput- 
ing the set of all the feasible positions (the orientations are fixed) of the objects, where certain 
constraints on their relative position are specified. The feasible region of (placements of) an 
object is called the free C-space of the object. The prohibited configurations of an object form 
a forbidden region. The C-space mapping for a particular spatial planning problem consists 
in partitioning the C-space into free and forbidden regions, where the latter are called C-space 
obstacles. See [29] and [63] for more details. 

*A full version of this paper is available at [42]. 

^PG-Consulting, Holon 58371, Israel. Email addreses: genap@post.tau.ac.il, 

gennadypus@walla . com. 



1 



1.1 Previous and related work 

Detailed surveys of previous work on spatial planning can be found in [9], [14], [21], [27], [29], 
and [63]. See also [2], [4], [26], [33], [54], and [58] and references therein for other related 
works on placement and spatial planning. 

The basic constraint on the relative position of geometric objects, used in spatial planning 
for generating the C-space maps, is the relative non-overlapping of objects; here the minimum 
translational distance between the objects must be greater than zero, or more generally, than 
some positive value. Indeed, this is a key requirement in robotics, packaging, and nesting. See, 
[4], [26], [27], [53], and [63]. However, the geometric problems arising in design and 
manufacturing require a placement of objects with more complex constraints on their relative 
position, such as constraints on the minimum and/or maximum translational distances between 
objects, and/or their Hausdorff distances. The problem of generating the C-maps with such 
complex constraints is also interesting theoretically. 

Placement problems taking into account the minimum translational distance (MTD) be- 
tween objects arise in industrial applications, concerning the cutting of materials, the layout of 
templates on a stock material, and the layout of an IC chip with geometric design constraints. 
See, e.g., [9], [29], [33], [53], and [54]. In these problems the MTD has to be at least the cutting 
tolerance of the machine that cuts the shapes out of stock material, or the minimal feasible dis- 
tances between electronic modules of an IC chip. Placement problems with constraints on the 
MTD between objects have been formulated in [53] and [54]. The papers [37] and [54] have 
considered placement problems with consrtaints on the minimal value of MTD between objects, 
and placement problems with constraints on both the minimal and maximal admissible values 
of MTD have been considered in [54]. Algorithms for solving various placement problems with 
constraints on the MTD have been considered in [37], [54], and [56] - [58]. 

The work in [39] has solved the problem of placement of a pair of objects with constraints 
on the minimum and maximum translational distances between them, and has also considered 
distances involving containment of objects. Placement problems with constraints on several 
types of translational distances between objects, including directional Hausdorff distances be- 
tween objects, have been studied in [40] and [41]. These works have also studied a Boolean 
function, called the geometric situation, which describes the system of constraints on the rela- 
tive position of objects, and have formulated and solved placement problems with several other 
types of geometric situations, namely, with rotational and dynamic geometric situations. 

An algorithm for computing the minimum Hausdorff distance between two planar objects 
under translation is given in [1]. (Efficient computation of Hausdorff distances has applications, 
e.g., in pattern recognition and computer vision; see [1] and [49].) This work has also studied 
placement problems taking into account bidirectional Hausdorff distances between objects. 

This paper presents a technique, based on the Minkowski operations, for generating the 
translational C-space maps for spatial planning problems with more general and more complex 
constraints on the relative position of geometric objects. This technique is an extention of that 
reported in [40] and [41]. The developed technique can also be used, respectively, for spatial 
planning with constraints on translational distances in a given direction, and/or on rotational 
distances between geometric objects, as well as for spatial planning with given dynamic geo- 
metric situation of moving objects. A full version of this paper is available at [42]. 

To formulate the problem let us first present the needed notations and definitions, and then 
consider the various standard distances between geometric objects. We assume that the ge- 
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ometric objects are regular sets (r-sets) in the Euclidean space _R", for n = 2 or 3, i.e., 
bounded, closed, and semi-analytic subsets of R" [43]. This means that, for any r-set A of 
i?", A = ki(y4), where i and k denote the interior and the closure of sets, respectively [43]. 
The complement and the boundary of A are denoted by A'^ and dA, respectively, and a copy 
of A translated by a point (or a vector) p is denoted hy A + p. We denote by A^ a copy of A 
rotated by an angle 6 about the origin O [24]. The regularized set operations on two objects A 
and B in i?" are defined A®* B = ki(y4 ® B), where ® e {U, n, \}; see [43] and [44] for 
details. The regularized complement A'^* of A is defined as A'^* = ki(A'^). The r-sets are not 
algebraically closed under the standard set operations, but they are closed under the regularized 
set operations [43]. For example, the standard intersection ACi B of r-sets A and B needs not 
be regular, since its boundary may have dangling faces/edges and/or isolated points, but An* B 
is an r-set. 



1.2 Standard distances between geometric objects 

Let us consider the following distances between objects A and B (see [1], [19], and [27]): 

di{B, A) = inf inf \\a — b\\; d2{B, A) = sup sup \\a — b\\ ; 

a€AbeB aeA b€B 

h{A, B) = sup inf \\a — b\\; h{B, A) = sup inf \\a — b\\; 
aeA beB Q^eA 

d*{A,B)= inf^inf ||a* -fell; H {A, B) = max{h{A, B) , h{B , A)} , 

where || ■ || is the Euclidean norm. The distance h{A, B) is called the directed Hausdorff distance 
from A to B, and the distance H{A, B) is called the Hausdorff distance between A and B, 
respectively; see [1]. H{A, B) is a metric on _R". 

These distances between r-sets have the following basic properties [19]: 

di{A, B)>0 AnB = d2{A, A) = diam(A); H{A, B) = A = B. 

Clearly, d*{A, B) = diiA", B). Then we have d*{A, B) > BcA. 



1.3 Problem formulation 

Let A be an unmovable object, and let B be another object, allowed only to translate. Let us 
consider the following problem: 

Problem I For given Ai, . . . , Ae, and the corresponding constraint 

ui{B + p,A) = |[rfi(fi + p,A) < Ai] 01 [d,{B + p, A) > \2\] 
©2 { + P,A)< As] 01 [d^iB + p,A)> A4] } 

02 { [d*{B +p,A)< As] 01 [d*iB +p,A)> Ae] }, 

where 0i(2) e {V, A}, find the region Ni{B, A) of all the feasible translations B + p of B,m 
which i'j{B + p, A) holds. 

Thus, our goal is to find all the feasible positions p of B with respect to A, under the above 
constraint on their relative position. This problem is generalization of the well known Findspace 
problem, formulated in [29]. If we let p = O, then the function vj^B, A) can be interpreted as 
the generalized Boolean distance query; see [27] for detailes. 
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2 Preliminaries 

2.1 Minkowski operations 

The Minkowski sum, and the Minkowski difference of objects A and B are defined as 

A®B = {a + b\ ae A,be B} = [j{A + b), and AeB = p| (A + 6) = {A^eBf, 

beB beB 

respectively [30], [49]. See Figure [2 a). The Minkowski subtraction is a dual operation of the 
Minkowski addition. Note that AQ)B = {A'^QBy. 




Figure 1 : (a) The Minkowski sum A(BB, and the Minkowski difference AQ B of objects A and 
B. (b) The objects A®B, and AQ B. Here pi G d{A®B), pa G d{iA®iB), {B + pi^2)(^A, 
and p3 G d{A fi), (i? + p3)cA, respectively. (Dashed lines show an objects B + p, where 
(a) p G OA, (b) p G d{iA®iB) (resp., G d{A Q B)). Dotted lines show pieces of d{iA®\B) 
(resp., d{\A®iB)).) 

Let 5 be the reflection of S with respect to the origin O, i.e., B = {—b \ b G B}. (For 
notational convenience, the object B is sometimes denoted by —B.) Then the dilation, and the 
erosion of A by B are defined as 

A®B = {a-b\ae A,be B} =[j{A-b), and AqB = f^{A - b) = {A'^By, 

beB beB 

respectively. See Figure (Ub). 

Many properties of the Minkowski operations are well known and well studied. See, e.g., 
[3], [14], [19], [26], [30], [31], and [49] for details. In this section we consider the properties of 
the Minkowski operations that we need for our purpose. 

In [8] and [14] it is shown that if an object A is a convex then AqB = AqCH{B) = 
AQext{B), where CH{B) denote the convex hull of B, and ext{B) denote the set of extreme 
points of B, i.e., the set of vertices of CH{B). In [31] it is shown that the Minkowski sum A®B 
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of two r-sets A and B always results in an r-set, whereas the Minkowski difference AqB could 
be a non-regular set; see Figure [IJa). 

Since A®{-{B +p)) = {A®B) -p, we have (see, e.g., [3], [17], [29], [30], [38], and [55]): 

{B + p)nA^^ p e A®B; 

(S+p)nA = ^ peiA^BY; (1) 
{B + p)r]A <^ p e d{iA®iB), 

where BnA = [{iA n 15 = 0) A {OA ndB 0)] denotes the outer touching of the objects 
A and B. The object A(BB is also called the C-space obstacle of B relative to A [29]. The set 
d{iA®iB) has been introduced in [3], where it is referred to as the outer envelope of A and B. 
In [3] it is shown that d{A®B)<Zd{iA®\B), and that the set d{iA®iB) may have coincident 
faces/edges and/or isolated vertices, which are removed from the open point set i{A®B)\ see 
Figured] 

From the relationships 9(Ae5) = diA^QiB) and Aq{-{B +p)) = (y4e5) -p, it follows 
that if AQB ^ 0, we get (see [30] and [39]): 

{B + p)cA ^^^peAeB; 

{B+p)(;tA ^ pe (AeBY; (2) 
{B+p)cA ^ pe d{AeB), 

where BcA = [{A" n 15 = 0) A {OA r\dB ^ 0)] denotes the inner touching of A and B. See 
Figure Kb). 

By observations of [3] and [31] we have {\A®\B)C\{A®B); {iA®\B)Z)\{A®B) and 

d{A®B)(Zd{iA®iB)- d{A®B) = d{iA®iB) = (9(A'=©i5); 

A®B = {iA®iB) U d{iA®iBy, A®B = i{A®B) U d{A®B)- (3) 
{A®BY = {iA®iBY\d{iA®iB); {A®BY = k[{A®BY]\d{A®B). (4) 
From the observations of [19] it follows that, for r-sets, we have 

A®iB = iA®B = iA®iB; A®B = k{iA®iB); A®B = A®iB = iA®iB. (5) 

(Clearly, for non-regular point sets, the above equalities do not necessarily hold.) 

Let us consider the difference B®A. From the properties of the Minkowski difference 
(see [30]) it follows that p e B®A if and only if p ^ A®B'', and then A n (5^= + p) = 0, i.e., 
Ac{B + p). Thus, B + p covers A if and only if p G B®A. In other words, 

Ac{B + p) ^r^peB®A; 

A^{B+p) ^ pe {B®AY; (6) 
AciB + p) ^ ped{B®A). 

(Note that, by previous observations, B®A, in general, is a non-regular set.) 

In case where both A and B are allowed to translate, we have {A ± q)®{—(B ± pY = 
{A®B) T p ± g, (A ± g)e(-(5 ± p)) = {A®B) T P ± q, and {-{B ± p))e(A ± q) = 
{B®A) ^p±q, respectively. Then all the relationships of ([T]), ©, and Q can be reformulated 
to handle this more general form. For instance, (B + p)c{A + q) <^=^ (p — g) G A®B. (For 
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A®B this well known fact (see, e.g., [17] and [49]) has widely been used in [2], [9], and [33], 
for solving various containment problems.) Clearly, for p = 0,we have 

BnA^ds ^ OeA®B; 
BcA ^ O G AeB; 

AcB ^ O G BqA, 

which are alternative formulations of the overlapping of the objects A and B, of the containment 
of B in A, and of the covering of A by i?, respectively. 



2.2 Distances between geometric objects concerning their outer relative position 

The standard minimum distance di{B, A) does not take into account the "amount" of intersec- 
tion between objects A and B, since di{B, A) = 0, for A D B ij}, regardless of how much 
they overlap. Minimum translational distance constraints that take into account penetration be- 
tween objects have been proposed in [5], [7], [35], [37], and [54]. The work in [42] consider 
one specific set of definitions of such minimum distances, since it has been defined in different 
ways. In this section we consider the translational distances, as defined in [7] and [37]. 
The minimum translational distance MTD(A, B), introduced in [7], is defined as 



MTD{A,B) 



-MTD+{A,B), for An 7^ 
MTD+{A, B), otherwise, 



where MTD+( A, 5) = inf{||t|| | An{B + t)}. 

The distances 71,2 A) between A and B, suggested in [37], are defined as 

f infceO(A®B) forAn5 = 0; 

71 (5, A) = < 0, forAnB; 

[ -infcG9(AeB) forAnB^i/}, 

72(5, A) = sup ||c||. 

ced{A(BB) 

See Figure m The distance 71 (-B, A) is defined by the above relationship only in case where 
d{AQ)B) = d{iA(BiB) (see subsection 2.1). Then, in general, we obtain that 



7i(S,A) = ( ^.^f;6a(iAeiB) 



inf^cflr;4m;R^ llc||, fov A H B = 

ell, otherwise. 



See Figure[3l By the observations of subsection 2.1, we get 7i(-B, A) = 7i(0, iA(BiB), and 

r < 0, for O G iA®iB; 
-fi{B,A)l =0, for O e d{iA(SiB); 
[ >0, forO G {A®BY. 

Note that in the above relationships the set iA (resp., iB) can be replaced by A (resp., B), since 
iA®iB = A®iB = iA®B, for r-sets. See [42] for more details. 

It can easily be shown that 7i(-B, A) = MTD{A, B). Therefore we denote the minimum 
translational distance between A and B by 71 (5, A). 
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(a) (b) (c) 



Figure 2: The distances 71 2(5, A), for various relative positions of A and B, in case where 
d{A®B) = d{iA®iB). Here (a) A n 5 7^ 0, (b) AnB, and (c) A n 5 = 0, respectively. 

7 1 (B, A) > 




V / 



(a) (b) W 

Figure 3: The distance ^i{B,A), for various relative positions of A and B, in case where 
c}(74©5)c9(iy4©ifi). Here (a) A fl i? 7^ 0; p^^, are the values of minimal translations of B 
in directions x and y, corresponding to 7i(-B, A), (b) Afi-B, and (c) A fl i? = 0, respectively. 



The properties of translational distances have been well studied. (See, e.g., [7], [17], [22], 
[34], [35], [37], [39], [54], and [56]. See also [42] for basic properties of 71,2(5, A).) In [7], 
[17], [37], and [56] it is shown that 'ji{B,A) (resp., 72(-B,A)) corresponds to the minimal 
(resp., maximal) translation B +pof B relative to A that reaches an outer touching {B +p)r\A, 
and that 7i(fi, A) = di{B, A), for A n 5 = 0, and 72(5, A) = 72(0, A®B) = ^2(5, A). The 
distances 7i,2(-B, A) are invariant with respect to both rotations and translations, i.e., 7i,2(-B^ + 
p, A'+p)= 71,2(5, A); see [35], [39], [56], and [57]. 

The papers [37] and [56] have considered the family of surfaces d{A(BB(B\K), for A > 0, 
where XK is the ball of radius A centered at O. In these works it is shown that 71 {B +p, A) = A, 
for p G d{A(BB(BXK). The surfaces d{A(BBQ\X\K), with similar properties, for negative 
values of A, have been defined in [56]. (Note that the above relationship holds only in case 
where d{A®B®XK) = d{iA®iBQ)iXK), for A > 0, and d{A®B) = d{iA®iB), otherwise.) 
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3 Parametric families of a single object and distances between a point and 
an object 

The parametric family of objects 



A®XK, for A > 0; 

AelXlK, for -TA < A < 0, 



where is the radius of the largest inscribed ball in A, is called the full parallel pencil of the 
object A [19], or, for A > 0, the offsets of A [47]. See Figuregta). The object Ti{-rA, K, A) is 
the locus of the centers of all largest inscribed balls in A, and, in general, is a (curved) face/edge. 
Clearly, ri(-rA,i^, A)cA. 

A second parametric family of objects 

r2(A, K, A) = XKeA, for A > Ra, 

where Ra is the radius of the smallest circumscribed ball of A, has been introduced in [40] 
and [41]. See Figure Hfb). From the definition of Minkowski operations, it follows that 
ri(A, K, A) = if and only if A < -r^, and r2(A, K, A) = if and only if A < Ra. Since 
XK is convex, it follows that r2(A, Ji, A) is convex, for any bounded A, and r2(A, K, A) = 
T2[X, K, CH{A)] = T2[X, K, ext{A)]; see subsection 2.1. The object T2{Ra, K, A) is a single- 
ton point, which is the center of the (unique) smallest enclosing ball of A. In general, the point 
of T2{Ra, K, A) needs not be in A; see Figure Sfb). From the definitions of ri(A, K, A) and 
r2(A, K, A) it follows that r2(A, K, y4)cri(A, K, A), for A > Ra- (An equality holds in case 
where A is a singleton point.) 



ri{\,K,A) 



■yi{p,A) 






(a) 



(b) 



Figure 4: The parametric families of objects Fi 2(A, A), for various values of A, and the 
distances 7i,2(p, A). Here r2(-RA,-^, ^) ^ A, dashed curves show dTi^2{X,K, A), and the 
dotted line shows a piece of dTi{X,K,\A). The objects Ti{—rA,K,A) and T2{Ra,K, A) 
consist of two singleton points, and a singleton point, respectively. 
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Let us consider the translational distances 71,2 (p, A) between a point p and an object A: 



miaedA\\p - a\\, for p e A"; 
7i(p,A) = { 0, forp e dA; 

-MaedA\\p- a\\, forpeiA, 

72 (p, A) = sup \\p - a\\. 

aedA 

To lack of space, all proofs here and in the following sections are omitted. The proofs of the 
observations, lemmas, and theorems are presented in [42]. 
The distances 71,2(^5 A) have the following properties: 

Observation 1 (a) For A > we have 

<X,forpeTi{X,K/iA); ( < X, for p e Ti{X, K, A); 

^i{p,A){ =X,forpe^T^{X,K,iA)■ 7i(P,^)< X, for p e [^T^{X, K,iA)] 
>X,forpe[T,{X,K,A)Y, I > X, for p e [T^{X, K,iA)] 



c. 

) 

c 



(b) For —va < a < we have 

( <X,forpeiT^{X,K,Ay, ( < X, for p e T,{X, K, A); 

-flip, A) I =X,forpedT,{X,K,Ay, 7i(p,^)< ^ X, for p e [dT^iX, K, A)Y; 

[ >X,forpe[T,{X,K,A)Y, [ > X, for p e k[T,{X, K, A)^. 

Observation 2 For A > i?^ we have 

( <X,forpeiT2{X,K,Ay, ( < X, for p e T^iX, K, Ay, 

^2{p,A) \ =X,forpedT2{X,K,Ay ^2{p,A) \ ^ X, for p e [dV^iX, K, A)Y; 

[ >X,forpe[T2{X,K,A)Y, [ > X, for p e kiT^iX, K, A)Y. 

It is clear that infpg/{n{7i(2)(p, A)} = —rA{RA)- Generally, for i = 1,2, we get 

rip A) I <A, forpGr,(A,ir,A); 

The topological properties of families Ti 2(A, K, A) have been studied in [42]. 



4 Correspondence between distances and the parametric families 

Let us consider the translational distance uj(B,A) between the objects A and B. We say 
that the distance u{B,A) corresponds to the parametric family of objects ^1{X, K, B, A) if 
and only if uj{B + p,A) < A, for p e Q{X,K,B,A). (Clearly, lj{B + p,A) > A, for 
p e [n{X,K,B,A)Y, and u{B, A) < A, for O G n{X,K,B,A).) The correspondence be- 
tween Lu{B, A) and fi(A, K, B, A) is denoted by iu{B, A) ~ fi(A, K, B, A). 

In special case where an object B is the origin O, the distance function cu{B + p. A) is 
reduced to the distance uj{0 + p. A) = uj{p. A) between a point p and an object A, and the 
family of objects VL{X, K, B, A) is reduced to VL{X, K, O, A) = Q{X, K, A). 

Consider the families of objects Ti 2(A, K, A), and the distances 71,2(0 + p,A). 
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Lemma 3 Fori = 1,2 we have: (a)ji{0,A) ~ Ti{\,K,A). (b)'yi{±p,A) =7^(0, A =Fp). 

Follow from the relationship (|7]), and since A(B{±p} = A^ p. 

Since diam(A) = 2Ra, it can easily be shown that 72(^5 A) = diam({j9} Li A) > A if 
and only if p e [r2(A, K, A)Y, for A > 2Ra. Then we get diam({0} U A) ~ r2(A, K, A), for 
A > 2Ra. (Note that 72 (p, A) < diam({p}UA) = diam(A) = A if and only if p E T2{\,K,A), 
for A < 2Ra.) 

The properties of the distances 7i,2(p, A) and their corresponding families Ti 2(A, K, A) in 
case where A = IJj=i and/or A = fXj^i Aj, for A 7^ 0, have been studied in [42]. 
In general case where 5 is a geometric object, but not a single point, we have 

Lemma 4 (a) ^i{B +p,A) = -fi{B, A-p) = -fi{p, iA®iB), fori = 1,2. (b) -fi{p, iA®iB) = 
'yi[B + a-p,A-{l-a)-p],forO<a<l. (c) ^i{B' ± p, A' ± q) = ^,{B' t q, A' t p) = 
-fi[±p T q, {iAmBYl for I = 1,2. 

Remark 1 For i = 2, the sets iA and/or iB can be replaced by A and/or B, respec- 
tively, since sup{||c|| | c E d{iA®iB)} = sup{||c|| | c E d{A®B)}, i.e., 72(p, iA©i5) = 
72(p,A©5). 

Remark 2 By the first of relationships ©, we have ji{B, A) = ji{B*, A*),fori = 1,2, 
where A* E {A, iA} and B* E {B, iB}, respectively. 
Let us first consider the parametric family of objects 



T^{X,K,B,A) 



{A®B)®XK, forA>0; 
{A®B)e\\\K, for - r^^s < A < 0, 



(see Figure [Sta)). Since ri(A, K, B, A) = Ti{X, K, A®B), the object Ti{-rj^^^, K, B, A) is 
the locus of the centers of all largest inscribed balls in A®B. 



Observation 5 (a) For A > we have 



<\,forpET,{\,K,B,A)- 
< \,forpETi{\,K,B,iA); 
= \, forp EOT i{\,K,B,iAy, 
>X,forpE[Ti{X,K,B,iA)r. 



< A, forp E Ti{X,iK, iB,iA); 

< A, forp E ill (A, iK, iB, iA); 
= A, forp E dTi{X, iK, iB, iA); 
> X, forp E k[Ti{X, iK, iB, iA)]". 

See Figures Oa) and [It a). Hence, the set 

is the surface of the function 71 (-B + p,A), and infpg/jn{7i(i? + p. A)} = — ri^^i^j. Thus, we 
get 




10 



Theorem 6 71 (5, A) 



ri(A,ir,B,A),/orA>0; 
ri(A, iK, iB, iA), for - r^^^i^ < A < 0. 

Note that in case where A > the Theorem [6] has also been proved in [37] and [56]. 

Denote by dp{A, B) the penetration depth of A and B (see [35]). In [42] it is shown that 

7i(5, A) = -dpiiA.iB), for A n ^ 0. Then we have dp{A, B) ~ Ti{X,K,B,A), for 



-rA(BB < A < 0. 



Let us next consider the parametric family of objects 

r2(A, K, B, A) = XKe{A®B), for A > i?^^^, 

proposed in [39] and [40]; see Figures Ob) andIMb). Since r2(A, K, B, A) = r2(A, K, A®B), 
the object T2{Raqb, K, B, A) is a singleton point, which is the center of the (unique) smallest 
enclosing ball of A(BB. The object r2(A, K, B, A) is convex, for any bounded A and B, and 

r2(A, K, B, A) = r2[A, K, CH{A®B)] = T^iX, K, extiAoB)]. 



ri(Ai, if, B, A) 



t 



S(A©S) 



B + p 



7l(B +p, A) 




P, A) 



Figure 5: The objects ri,2(Ai,2, K, B, A). Here 71,2(5 + p,A) = Ai,2 and Ai > 0. A dotted 
line shows a piece of dTi{Xi, K, B, iA), and 7i(-B + q,A) = Ai, for q E dTi{Xi, K, B, iA). 



Observation 7 For X > R^^p we have 



72(5 +p, A) 



<XJorpET2{\K,B,A)- 
< X,forpeir2iX,K,B,A); 
= X,forpedT2{X,K,B,A); 
>X,forpek[T2{X,K,B,A)Y. 



(Clearly, infpgRn{72(i? + p, A)} = Ra^b-) Then we get the following: 

Theorem 8 72(5, A) ~ r2(A, K, B, A), for X > Ra^^b- 

By previous observations, we also obtain that 72 (-B + p, A) = diam((i? + p) U A) > Xif 
and only ifp e [r2(A, K, B, A)Y, for A > 2/2^^^. Hence, diam(5 U A) ~ r2(A, K, B, A), for 



A > 2Ra(qb- 



The additional (topological and set theoretic) properties of the distances 71,2(5, A) and the 
families ri,2(A, K, B, A) have been studied in [42]. 
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Figure 6: The objects (a) ri(Ai, iK, iB, iA) and (b) r2(A2, K, B, A). Here 7i,2(S+p, A) = Ai,2 
and Ai < 0. A dashed lines show d{iA(BiB). 

5 Distances between geometric objects concerning their inner relative posi- 
tion 

The duality of the Minkowski operations provide estimating the inner relative position of geo- 
metric objects. 

5.1 Distances concerning the containment of objects 

Let us consider the following translational distances, introduced in [39] and [40]: 



See Figure|7l (The distances ?7i,2(-B, A) are defined only in case where AqB ^ 0.) 

The properties of the distances ?7i,2(5, A) have been studied in [39] and [41]. The distance 
r]i{B, A) (resp., 'r]2{B, A)) corresponds to the minimal (resp., maximal) translation B + p of B 
relative to A that reaches an inner touching {B +p)(zA. Since d{AQB) = d{A^(BiB), we have 
7]i{B, A) = -7i(B, A"") and r]i{B, A) = -d*{A, B), for BdA. 

Lemma 9 (a)T]i{p,A) = -^^{p, A),fort = 1,2. (b)rji{B+p,A) =r]i{B,A-p) =7]i{p,AQB). 
(c) rii[B + a ■ p,A - (1 - a) ■ p] = r]i{p, AqB), /or < a < 1. (d) rii{B^ ±p,A^ ±q) = 
ri.iB'^ Tq.A^ Tp) =V^[±PTq,iAeBy],for^ = 1,2. 

It can easily be shown that 




for B A; 
for Be A; 
for Be A, 



V2{B,A) 



sup 

cea(AeB) 



c 




<0, forOGi(Ae5); 
= 0, for O e diAeB); 
> 0, for O e {AQBy. 
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))i(B,A)<0 





AQB 





1)2(0, A e 5) 



Vi(0,AeB) 



(b) 



Figure 7: The distances 771 2(-B, A), for various relative positions of A and B. Here (a) BcA, 
(b) -BcA, and (c) 5^/1, respectively. 



See Figure|71 The distances 771 2 A) have used in [40] and [41] to describe the constraints on 
the relative position of objects in containment problems. 

Let us consider the following parametric families of objects [39], [40]: 



H,{\,K,B,A) 
H2{X,K,B,A)-- 



{AeB)®XK, forA>0; 
{AqB)q\X\K, for - r^QS < A < 0, 

\Ke{AeB), forA>i?^e^, 



See Figure [H 



d{AQB) 



B + p 



Hi(Xi,K,B,A) 



7ji(B +p,A) 



(a) 




Figure 8: The objects -ffi, 2(^1,2, B, A). Here 77i,2(-B + p, A) = Ai,2 and Ai < 0. 
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Observation 10 (a) For X> Owe have 



r]i{B+p,A) 



< A, forp E Hi{\,K,B,A); 

< \,forp e Hi{\/iK,B,A); 
= \,forpedHi{\,iK,B,A)- 
>\,forpE[H,{\,iK,B,A)Y. 



r]i{B + p,A) 



(b) Fori = 1, 2, —r^Q^ < Ai < 0, and A2 > RaqB' respectively, we have 

< XiJorpEHi{\,K,B,A)- 

< Xi,forpe iHi{Xi,K,B,A); 
= Xi,forpedHi{Xi,K,B,A); 
>Xi.,forpek[Hi{X„K,B,A)Y. 

Then mfpg^n {771(2) (-B + p, A)} = —ry^QBi^AeB)- Thus, by above, we obtain the following: 
Theorem 11 r],{B, A) ~ H,{X, K, B, A), for i = l,2. 

The additional properties of the distances 7712(5, A) and the families Hi 2{X, K, B, A) have 
been studied in [42]. 



5.2 Distances concerning the covering of objects 

In this section we introduce the distances 5i,2(-B, A) involving the covering of A by B: 



5iiB,A) 
52{B,A) 



infce9(BeA) l|c||, for A ^ fi; 



0, 



fovAcB: 



-infcGa{BeA) l|c||, fovAcB, 
sup ||c||. 



See FigurelU (The distances 5i^2{B, A) are defined only in case where BqA ^ 0.) 

<5i(B,A) = o y 





BQA 





71(0, Be A) 



Si{B,A) > 



Figure 9: The distances 61^2(8, A), for various relative positions of A and B. Here (a) AcB, 
(b) AgB, and (c) A <;t B, respectively. 



The properties of the distances 5i,2(-B, A) and 77i,2(-B, A) are similar. The distance Si{B, A) 
(resp., 52(5, A)) corresponds to the minimal (resp., maximal) translation B + p of B relative 
to A that reaches an inner touching Ac{B + p). The following helpful properties of (5i,2(5, A) 
hold: 5i{B, A) = -71 (5^ A) and 5i{B, A) = 7]i{A, B), for z = 1, 2. 
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Lemma 12 (a) Si(B + p,A) = Si{B, A - p) = Si{BeA, -p) = -fi{p, BeA), fori = 1,2. 
(b)5i[B + a-p,A- {l-a)-p]= SiiBeA, -p), for < a < 1. (c) 5^{B^ ± p, A^ ± q) = 
S,{B' Tq,A'Tp) = S,[{BeAY, TP±q]= 1^[±P T q, {BQAY], for z = 1, 2. 

It can easily be shown that 

r < 0, forO G i{BeA); 

6i{B, A) I =0, for O e d{BeA); 

[ >0, for O G {BeAy. 

See Figure|9l The distances 5i,2(5, A) can be used to formalize the constraints on the relative 
position of objects in covering problems. 

Let us consider the parametric families of objects 



A^{X,K,B,A) 



{BeA)®XK, forA>0; 
{BeA)e\X\K, for - r^QA < A < 0, 



A2iX, K,B, A) = XKe{BeA), for X>R^QA. 
See FigurelH Since A,(A, K, B, A) = -Hi{X, K, A, B), for i = 1, 2, we get 

Observation 13 (a) For A > we have 

( <X,forpeA,{X,K,B,Ay, 
;s(RMr.A\} < Kforpe A^{X,iK,B,Ay, 
5^{B + p,A)^ =X,forpedA,iX,iK,B,A); 

[ >X,forpe[A^iX,iK,B,A)Y. 



(b) Fori = 1, 2, — r^g^ < Ai < 0, and X2 > RbqA' respectively, we have 

( < Xi, for 

6,{B+p,A) 



< X„forpeA,{Xi,K,B,Ay, 

< Xi, forp eiAi{Xi,K,B,A); 
= Xi,forpedA,{Xi,K,B,A); 
>\,forpek[\{\,K,B,A)Y. 



Thus, infpg/jn{5i(2)(5 + p, A)} = — ?"BeA(-^BeA)- Hence, we can conclude 
Theorem 14 6i{B, A) ~ Ai(A, K, B, A), for i = l,2. 

Note that the Theorem [14] also follows from the claim (a) of Lemma [I2l and from the relation- 
ship Ai(2) (A, K, B, A) = ri(2)(A, K, BeA). 

The additional properties of the distances 5i,2(-B, A) and the families Ai^2(A, K, B, A) have 
been studied in [42] . 



6 Hausdorff distances and corresponding families of objects 

Since H{A,B) = inf{A > | BcTi{X, K, A), AcTi{X, K, B)} [19], we have A) = 
inf{A > I BcTi{X,K,A)}. Then h{B,A) = in case where BcA, i.e., the distance 
h{B, A) does not take into account the "amount" of containment of B in A. 
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Figure 10: The objects (a) Ai,2(Ai,2, K, B, A). Here 5i,2(5 + p,A) = Ai,2 and Ai < 0. 



In [41] has been introduced the signed distance ^{B, A) = supfeg^{7i(6, A)}, eliminating 
this shortcoming, and it is shown that 

„rR A^-S ^^B,A), fori? ^ A; 

\ Vi{B,A), otherwise. 

Then the signed distance fi{A, B) can be defined as 

' ^ \ 8i{B,A), otherwise. 

Hence, H{A, B) = max{/i(A, B), ^{B, A)}. 

In [1] and [41] has been suggested the parametric family of objects Mi{\, K, B, A) = 
Ti{\,A)eB, for A > -r^^^, and it is shown that ij.{B,A) ~ Mi{\, K, B, A); see Fig- 
ure [nja). The family of objects M2(A, K, B, A) = Ti{X, B)eA, for A > 0, has been intro- 
duced in [1], and it is shown (in our notation) that h{A, B) ~ M2(A, K, B, A), and H{A, B) ~ 
M3(A, A', 5, A) = [Mi(A, K, B, A) nM2(A, K, B, A)], respectively. See FigureiHb). Clearly, 
for -r^QA < A < 0, we have ij.{A, B) ~ M2(A, B, A). 

For notational convenience, the distances ix{B,A), fi(A,B), and H{A,B) are sometimes 
denoted by mi{B, A), m2{B, A), and m^{B, A), respectively. Then we get 

Theorem 15 mi{B, A) ~ Mi(A, AT, 5, A), for i = 1 - 3. 

The objects Mi_3(A, K, B, A) have the following simple properties: 

Mi(2)(A,i^,5,A) = -A/2(i)(A,A;A,5); Ms{\,K,B,A) = -M,{X, K, A, B). 



Observation 16 (a) For i = 1 — 3, and A > we have 

( <\,forpeMi{\,K,B,A); 

m (B + vA)l < ^' /^^P ^ ^^); 

rn^yn ^ p, ^) ^ ^ ^ ^ ^^^^^^ A)\iMi{\ iK, iB, iA); 

[ >\forpe k[Mi(A, iAT, i5, iA)f. 
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Figure 11: The objects (a) Mi^2{^i,2, K, B, A). Here, respectively, fi{B + p,A) = Ai > 
0, fi{A,B+p) = A2 > 0^ and dTi{\i,K,A) = dTi{\i,iK,iA), dTi{\2,K,B) = 
dT,{X2,iK, iB). 



(b) For i = 1,2, — r^eij < Ai < 0, and — T^eA < A2 < we have 

( < X^,forpeMi{Xi,K,B,A); 
rr,(R^r^A\} < M P ^ iMi{Xi, K , B , A); 
^ ^' M = A„ forp e 9M,(A„ K, 5, A); 

>K,forpek[Mi{Xi,K,B,A)Y. 



3 In [42] it is shown that, in contrast to the distances considered in Sections 2 - 
where the distance n{B + p, A) (resp., n{A, B + p)) is equal to A may have the 



Remark 3 

5, the region 
non-empty interior. 

The additional properties of the Hausdorff distances and their corresponding families have 
been studied in [42]. 



7 Translational distances between geometric objects and translational geo- 
metric situations 

The distances considered in Sections 2 - 6 are referred to as the translational distances (or T- 
distances). Let TV{B,A) = {71,2(5, ^71,2(5, ^1,2(5, A), mi„3(5, A)} be a collection 
of the T-distances between B and A, and let u{B, A) G TV{B, A). The T-distances uj{B, A) 
and their corresponding families ^1{X, K, B, A) have the following properties: 

uj{B + p,A) = uj{B,A-p) = uj[B + a-p,A- (1 - a) ■ p], forO < a < 1; 

uj{B' ±p,A'±q) = uj{B' ^q,A'Tp); 
fi(Ai, K, B, A)Cfi(A2, K, B, A), for any bounded A and B, and Ai < A2; 
fi[A, K,B + a-p,A-{l-a)-p] = fi(A, K, B, A) - p, for < a < 1; 
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n{\, K, ±p,A^ ±q)= K, ^q^A^Tp) = [^(A, K, B, A)Y TP±q- 

From the last relationship it follows that uj{B ± p, A ± q) < A, for ±p T Q ^ K, B, A). 
By the previous observations, we get 

(A n 5 = 0) A (5 ^ A) ^ [71 (5, A) > 0] A I [viiB, A) > 0] V [^B, A) > O] |; 

(A n 5 ^ 0) A (AcB) ^ [71(5, A) < 0] A I Hb, A) < 0] V [KA 5) < 0] }. 

Hence, the various situations of the relative position of objects A and B can be described by the 
system of constraints on the T-distances between A and B. 

Definition 17 The relationship u{B + p, A) = [ijj{B + p,A) Q A], where G {<, =, >}, is 
called the primitive translational geometric situation (PTGS) of an object B with respect to an 
object A. 

Let us consider the PTGS u{B + p, A) as an event. Then the class of all the possible 
PTGS's S^{B, A) forms an algebra At of events [23]. For S^{B, A) permitting the following 
definitions: 

1. The union of events: i'i{B + p, A) V i'2{B + p, A). 

2. The intersection of events: i'i{B + p,A) A h'2{B + p,A). 

3. The complement of event: [z/(i? + p, Ayf. 

4. The certain event I is the union of all the PTGS's in S'^^B, A). 

5. The impossible event is an impossible relative position of objects, for given z/(i? +p, A) . 
The class S{B,A) of PTGS's, comprising S^{B,A) and 0, forms a completely additive 

Boolean algebra; see [23]. 

Let us consider the union of PTGS's. (Recall that u{B, A) e TV{B, A).) Then we get 

[uj{B + p,A) < A] = [lu{B + p,A) < A] V [lu{B + p,A) = A]; 

[u{B +p,A)j^\] = [u{B + p,A)<\]y [u{B + p,A)> \]; (8) 

[u{B + p,A) > A] = [uj{B + p,A) > A] V [u{B + p,A) = A]. 

(Clearly, [u{B + p, A) < A] A [u{B + p,A) > \] = 0.) We define also the TGS of type 
[uj{B + p, A) min], since in Sections 4-6 have been obtained the possible minimal values 
of the T-distances. We next add these TGS's to the set of PTGS's: 

Definition 18 The relationship v{B + p,A) = [lo{B + p,A) Q X or uj{B + p, A) min], 
where G {<,<,=, 7^, >, >}, is called the basic translational geometric situation (BTGS) of 
B relative to A. 

Since [uj{B + p,A) < A]^ = [lu{B + p,A) > A], [lu{B + p,A) = A]^ = [uj{B + p,A) ^ A], and 
[lu{B+p, A) > XY = [lu{B+p, A) < A], then, by the relationships of ([8]), and, by DeMorgan's 
laws, we obtain the following useful relationships: 

[uj{B + p,A) < A] = [io{B+p,A) < A] A [uj{B + p,A) ^ A]; 

[u{B +p,A) = X] = [iu{B + p,A)<X]A [lu{B + p, A) > X]; 
[uj{B + p,A) > A] = [uj{B+p,A) > A] A [uj{B + p,A) ^ A]. 
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Definition 19 The Boolean function u{B+p, A) = f[ui{B+p, A),..., Uk{B+p, A)] ofBTGS's 
i^i{B + p, A), for i = 1, . . . , k, is called the translational geometric situation (TGS) of B with 
respect to A. 

A particular TGS describes the constraints on the relative position of objects, and can be used 
as the objective function in the spatial planning problems. 

8 Constructing the feasible region of an object for given translational geo- 
metric situation 

Let f2(A, K, B, A) is the corresponding family of the distance uj{B, A). From the observations 
of Sections 4-6, and from the Definition [T8] it follows that, for given BTGS ^{B + p,A), all 
the feasible translations B + p of B with respect to A are obtained if and only if p belongs to 
the region N{B, A), where u^B + p, A) is true. For example, 

ui{B + p,A) = [u{B+p,A)<\i] ^ peNi{B,A) = n{Xi,K,B,Ay, 
U2{B + p,A) = [uj{B +p,A)>X2] ^ pe N2{B, A) = [^(Aa, K, B, A)Y. 

The correspondence between u{B +p, A) and N{B, A) is denoted by u{B +p, A) ~ N{B, A). 

Consider the TGS u{B+p,A) = i^i{B +p, A) Au2{B +p, A), where the BTGS' s vi^2{B + 
p, A) are as above. Clearly, its corresponding region is N{B, A) = Ni(B, A) n N2{B, A). In 
case where Ai < A2, we have N{B, A) = 0, and therefore u{B + p, A) is an impossible TGS, 
i.e., u{B +p,A) = 0. For the TGS u{B + p,A) = vi{B + p, A) V i'2{B + p, A) we have 
N{B, A) = iVi(5, A) U N2{B, A). If Ai > A2 then N{B, A) = i?". (Note that in this case 
iy{B+p,A)^I.) 

Consider next the TGS 1/(5 +p, A) = [lv{B+p,A) < Xi]y[uj{B+p, A) < A2], where Ai < 
A2. Then N{B, A) = n{Xi, K, B, A) U ^(Aa, K, B, A) and f](Ai, K, B, A)cn{X2, K, B, A), 
forAi < A2. Hence, N{B, A) = n{X2, K,B, A). Thatis, theTGS A) corresponding to 

the region N{B, A) can be represented in different ways, e.g., as iy{B+p, A) = [uj{B+p, A) < 
X2] or ^{B+p, A) = \/^_^[uj{B+p, A) < Aj], where Aj < A2, for z = 1, . . . , n. Thus, in general, 
the TGS ^{B + p, A) does not unique relative to its corresponding region N{B, A). However, 
for given N{B, A), the unique (minimal) corresponding TGS u{B + p, A) can be constructed. 
Finally, if in the above example, we let N{B, A) = Q{X, K, B, A) U Q{X2, K, B, A), where 
Ai < A < A2, then we obtain that the TGS [uj{B + p,A) < Ai] does not corresponds to the 
region ^1{X, K, B, A), however z/(i? + p, A) ~ N(B, A) in this case. Therefore, we get the 
following: 

Proposition 20 (a) The TGS u{B +p, A) = f[ui{B +p, A), . . . , u^iB +p, A)] corresponds to 
the region N{B, A) = F[Ni{B, A),..., N^iB, A)] if (but not only if) u,{B+p, A) ~ Ni{B, A), 
for i = 1, . . . ,k, and the operations of union, intersection, and complement of the BTGS 's 
UilB+p, A) ofv{B+p, A) correspond to the operations of union, intersection, and complement 
of regions Ni{B, A) ofN{B, A), (b) ^{B+p, A) is an impossible TGS if and only ifN{B, A) = 
0. (c) v{B + p, A) is a certain TGS only if (but not if) N{B, A) = i?". 

The region N{B,A) is called the solution of the TGS i/(i? + p,A). (Note that the region 
Ni{B, A) of N{B, A) is solution of the BTGS UiiB + p, A), for i = 1, . . . , k.) 
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Remark 4 Let ^{B + p, A) be the BTGS concerning the T-distance uj{B + p, A), and let 
uj{B +p, A) ~ K, B, A), and ^{B +p, A) ~ A^(i?, A), respectively. In the next paragraph 
we assume, for simplicity, that the region N{B, A) does not have coincident faces/edges and/or 
isolated points, which are removed from the interiors of the objects of r2(A, K, B, A), for any 
A. Then we get 



In this special case the set dQ{X, K, B, A) is the surface of the function uj{B + p,A). 

The types of objects satisfying the above assumption are sufficiently wide. For example, 
the convex objects, the polygons/polytops and/or curved objects in general position, i.e., which 
do not have parallel edges/faces. Thus, one can to construct the solution of a particular TGS 
according to the more simpler relationships than that is done in Sections 4-6. See [42] for 
details. 

Examples. Consider the TGS's h'i{B + p, A) and their solutions Ni{B, A), denoted by ui, Ni, 
for short. Denote also lj{B + p, A), Q{X, K, B, A) by u, Q{X). Below we let = 1,2. 



1 i^i = (7. < A.) A (7, > A,); N, = r,(A,) n k[r,(A,)]'= = r,(A,)\ir,(A,). 

la ly^a = < A,) A (7,- > A,); iVia = ir.(A,) n [r,-(A,)]' = ir,(Ai)\r,(A,-). 
lb z/ifc = (7^ < Xi) A (7,- > A,); iVi6 = r,(A,) n [r,(A,)]' = r,(A,)\r,(A,-). 
Ic z/ie = (7^ > A.) A (7, > A,); iVic = k[r,(A,)]'= n k[r,(A,)]^ = k[r,(A,) u r,(A,)]^ 

2 z/^ = (7* = A.) A (7, = A,); = dT,{X,) D dT,{X,). 
3 1^3 = (71 > Ai) A (72 ^ min), where Ai ^ r^^^^y, 
N, = k[T^{X^)r H T,iR^^^) = r2(i?^eB)\iri(Ai). 
4 1/4 = (?7i min); = Hi{-r^Q^). 
5 1^5 = (^1 < Ai < 0) A (r/2 ^ min); N5 = Hi{Xi) n H2{Raqb)- 
6 z/6 = [(71 = Ai) A (72 = A^in)] ^0; iVe = dT^{Xi) n dV^iX^,^), 



See Figure [I2l The analysis of the Examples 1-6 can be found in [42]. 

Remark 5 The region N{B, A) may have various topology. It can be open or closed, regular 
or non-regular, bounded or unbounded, connected or disconnected. So, in the above examples 
the region A^^i is closed bounded, whereas the region Nia is open bounded; the region Nic is 
closed unbounded. In case where N{B, A) contains the subset of its boundary, it is neither open 
nor closed. Thus, in general, N{B, A) is an object with non-manifold boundary. The region 
Nn, give an example of such an object. In case where i = j, Aj > Xj, and Xij > 0, respectively, 
the object A^i is regular, i.e.,A''i = Ti{Xi)\*Tj{Xj). (Note that in case where Aj = Xj we get 
Ni = dTi{Xi), i.e., the object Ni is non-regular.) 




where A 



'mm 



inf{A I r2(A)cri(Ai)}, for T^iRAeB) e ri(Ai); 
inf{A I r2(A)nri(Ai)}, forV^iRAs^B) i ri(Ai). 
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Figure 12: The corresponding regions Ni of TGS's z//, for / = 1,...,6, in case where 
V2{Ra®b) ^ ri(Ai). Here z = 2, j = 1, and Ai < A2, respectively. The region Ni is 
simply connected, N2 = {pi,P2}, and A^3 = T2{Raqb), respectively. The region A^4 is a curve, 

N5 = 0, and Nq = {gi, ga, Qs}- 



Let us next consider the following problem: 
Problem II Find the region Njj, corresponding to the TGS 

i^ii = [(Ai < 71 < A2) A (A3 < 72 < A4)] V [(A5 <Vi< Ae) A (A7 < r/2 < Ag)] . 
By Example 1, the solution is the region 

Nil = { [ri(A2)\iri(Ai)] n [r2(A4)\ir2(A3)] } u | [H,(x,)\iH,{x,)] n [H2{Xs)\iH2{x, 

See Figure [T3l Note that in case where Ai > A2 or A3 > A4, and A5 > Ae or Ay > As we have 
Nil = and un = 0. 

Let us turn to solve the problem, formulated in subsection 1.3. The following regions Nj, 
for ©1 = A, 02 = V, and N'/, for 0i = V, 02 = A, are the solutions of the Problem I: 



N'j = [ri(Ai)\iri(A2)] u [r2(A3)\ir2(A4)] u [H^{x,)\iH,{x,) 
N'{ = [ri(Ai) uk[ri(A2)]^] n [r2(A3) uk[r2(A4)]^] n [iJi(A5) uk[iJi(A6)]^]. 

Clearly, if Ai < A2, A3 < A4, and A5 < Ae, then Nj = ^ and = 0. In case where Ai > A2, 
A3 > A4, and A5 > Ae we have N'/ = R^. (However, v'l 7^ / in this case.) 
The additional properties of the TGS's have been studied in [42]. 



9 Applications 

In this section we consider the spatial planning problems with more general and more complex 
constraints on the distances between geometric objects. We also briefly consider the several 
other types of geometric situations: the translational geometric situation in a given direction, 
the rotational, and the dynamic geometric situations, respectively. 
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Figure 13: Illustration for the Problem 11. Here the region Nu is 5-connected. 



9.1 Findspace problem 

Let A = {Ai, . . . , An} be a collection of n, possibly intersecting, obstacles Ai, completely 
contained in the region R, and let B be the object moving relative to A under translations. See 
Figure [T4T a). (For notational convenience, we also denote A = ljr=i ^j-) 

The translational Findspace problem [29] is to find all the possible translations {B + p) cR, 
such that {B + p) f] A = $An this case p is called the free position. (If (B + p)r]A, then p is 
called the semi-free position [4].) 

The C-space obstacle of B relative to A is defined as COb{A) = {p \ {B + p) H A 
0} = A®B = \JUCOB{Ai) = \JU{^i®B) [29]. (The object [COb{A)Y is called the 
free C-space of B relative to A.) The C-space interior of B relative to the region R is defined 
as CIb{R) = {p I {B +p)CR} = [COb{R^)Y; see Figure USa). By the definition of the 
Minkowski difference, we get CIb{R) = RQB. 

In [14] it is shown that the set of all the feasible positions of B relative to A and R can 
be represented as FPb{A,R) = {R\A)eB = {ReB)\{A®B) = CIb{R)\COb{A)\ see 
Figure Oa). 

Let us consider the Findspace problem in terms of the translational geometric situations 
(TGS), studied in Sections 7 and 8. The conditions {B + p)cR and (B + p) f] A = $ can be 
formalized as [r]i{B + p, R) < 0] and [71 (-B + p,A) > 0], respectively. In [42] it is shown 
that [71 (-B + p, A) > 0] = /\i^i['JiiB + p, A^) > 0]. Then the Findspace problem can be 
reformulated as follows: Find the region N{B, A, R) corresponding to the TGS 

r^{B + p, A, R) = [r]^{B + p, i?) < o] /\ | /\ [71(5 + p, A,) > o] |. 

i=i ^ 
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Figure 14: The translational Findspace problem in case where A = {Ai, A2, A^,}. (a) Here 
the region FPb{A,R) is connected, and ]?i 2 are the free positions of B. Dashed lines 
show dlCOsiA)] and d[ClB{R)], respectively, (b) Illustration for the Problem III. Here 
Njjj{B,A, R) is a two-connected region and Ai = 0, A2 = |A/j|,and A3 = 1.5- A2, respectively. 
Dashed curves show (9[C0b(A^, A)] and ^[^/^(Aij, R)]. 



Since [71(5 + p,Ai) > 0] = [71(5 + p,A,) < Of, we have ALi bii^ + p,A,) > O] = 
{ Vr=ibi(-^ + — ^]Y- Hence, the solving of the Findspace problem can be represented 

as follows: 



N{B, A, R) = ifi(0, K, B,R)f]lf] [ri(0, A', 5, A 

^ 1=1 

n 

= H,iO, K,B,R)f][[j ri(0, 5, A,) 

1=1 

n 

= H,{0,K,B,R)[] T,(o,K,B,[jA 



= Hi{0,K,B,R)\T,{0,K,B,A). 

Clearly, i/i(0, K, B, R) = CIb{R) and ri(0, K, B, A) = COb{A). 
We next consider more general Findspace problems. 
Problem III Find the corresponding region Njjj{B, A, R) of the TGS 

/\[7i(fi + p,A,)>A, >0 
i=i 

Solving the Problem III. For given TGS i^ni{B + p, A, R), the C-space obstacle depends on 
= {^i}iLv corresponds to the TGS Vr=ibi(-^ + ^i) ^ -^j]' therefore it can be 
represented as COb(A^,^) = IJ"^^ ri(Ai, -B, Aj). The interior C-space corresponding to 
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the TGS [r]i{B+p, R) < \r < 0] is the region CIb{Xr, R) = Hi{\r, K, B, R). (See Sections 
4, 5, and Figure [T4lb).) Thus, we get 



Njn{B, A, R) = H^{\r, 5, i?) f| | f| [ri(A„ K, B, A. 

^ 1=1 



= H,{Xr, K, B, R)\ [ [j T,{Xi, K, B, A, 

1=1 

Problem IV Find the solution Njv{B, A) of the TGS 
Piv{B+p,A) = -fi,{B+p,A) < Ai] /\ ^-fi,{B+p,A) > X2 , where /i(2) = 1,2; Ai(2) > 0. 
Solving the Problem IV. The general solution is the region 

NjviB, A) = Ti^ (Ai, K, B, A) f| [r^, (A2, K, B, A)] ' = Ti, (Ai, K, B, ^)\r;,(A2, K, B, A). 
By observations of [42] and of Section 8, in case where /i(2) = 1 and Ai(2) > 0, we have 

ujviB + p,A) = \ \/[^iiB + p,A,)<\^'\ l/\ /\[^,iB+p,A,)>\ 
^ i=i ^ ^ i=i 

n n 

NiviB, A)=[[j ri(Ai, K, B, A,)] \ [ U ri(A2, K, B, A,)] , 

i=l i=l 

whereas, for /i(2) = 2 and Ai(2) > mmi<i<n{RAi®B}^ we get 

uiv{B + p,A) = l/\ [72(5 + P, A,) < Ai] I /\ I V [72(5 + p, A,) > A2] |; 

n n 

Niv{B, A)=[f] r2(Ai, K, B, A,)j \ [ f| r2(A2, K, B, A, 



i=l 



i=l 



See Figureini Note that uiv{B + p,A) = and Niv{B, A) = 0, for Ai < A2. 

Finally, we consider the Findspace problem concerning the covering of objects. 

Problem V Find the region Nv{B +p, A) of all the possible coverings of object A by object 
B + p, for given TGS 



uv{B+p,A) = l \/[52iB+p,Ai)<E,\ <5i{B+p,Ai)<\i <0 

^1=1 ^ ^ i=l 

The solution is 

Nv{B + p,A)= [[jA2{e,,K,B,A,)'\f]l f][A,{Xt , K, B, A,)\iA,{X- , K, B, A,)] I. 



i=l 



i=l 



For more simple constraint i^y(-B+p, ^) = /\^^^[5i{B +p, Ai) < Aj < 0], we get A'"y(_B, ^) 
nr=i Ai(Ai, K, B, A,i)- see Figure[ 
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Figure 15: Illustration for the Problem IV. Here A = {Ai, A2}, both Ai and A2 are line seg- 
ments, B = {O}, Ai > A2, and the region Njy{B, A) is connected. Dashed curves show pieces 
of arz,(i2)(Ai(2), 5, Ai(2)). Here (a) h(2) = 1, and (b) h(2) = 2. 

9.2 Placement of geometric objects 

An approches commonly used for solving the placement problems are the sequential-single 
method, suggested in [53] and [57], and the multiple placement of several objects, suggested 
in [2], [9], and [33]. 

The sequential-single placement consists in the sequential placement of the objects of A 
with respect to the container Aq in a fixed order, say Ai, ^2, . . . , A^, according to the special 
objective function. For instance, the valid position pi of Ai must be a point with extremal (or 
specified) values of coordinates. (In case of planar problems pi = (xj, yi) can have, for example, 
a minimal, maximal, or specified Xi and/or yi.) 

The multiple (or simultaneous) placement, as follows from its name, is independent on the 
order of placement of the objects of A relative to Aq, and provides the placement of each object 
of A, say Aj, taking into account the possibility of the placement of all another objects {Ai}, 
where 1 < i < n, and i ^ j. The goal is to find the set Pqj of all the feasible positions of Aj, 
for J = 1, . . . , n, with respect to Aq. 

The generalized sequential- single and multiple placement problems (i.e., the problems with 
more complex constraints on the relative position of objects than in [33] and [57]) and their 
solutions have been studied in [42]. 

9.3 Application to the other types of geometric situations 

The work in [41] has studied the following types of geometric situations: 

The translational geometric situation in direction u describes constraints on translational 
distances in a given direction between geometric objects. The minimum and maximum dis- 
tances 7i 2(-B, A, u) taking into account the outer position of B relative to A in direction u have 
been introduced in [22] and [37]. The minimum and maximum distances 771 2(-B, A, u) taking 
into account the inner position of B relative to A have been proposed in [39] and [41]. The 
parametric families of objects corresponding to the distances in a given direction are obtained 
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Figure 16: Illustration for the Problem V. (a) An objects A = {Ai, A2, A3} and B. (b) Here 
the region Ny{B, A) is siply connected, Ac{B + p), for p E Ny{B, ^); Ai = 0, A2 < 0, and 
A3 = A2, respectively. Dashed lines show 9Ai(Ai_3, 5, A). 

by using the partial vector operations on objects, which are generalizations of the Minkowski 
operations. See [41] and [45] for more detailes. 

The rotational geometric situation describes constraints on minimum and maximum rota- 
tional distances between geometric objects. Denote by A* the image of A in polar coordinates, 
i.e.. A* = {(r, 9) \ (r cos 9, r sin 9) E A}. Then a copy A'^ of A rotated by an angle (p around 
the origin O corresponds to a copy A* + {(0, (f))} of A* translated by a point (0, (p) in polar 
coordinates. Therefore one can use the partial vector operations in polar coordinates for mod- 
eling the relative position of geometric objects under rotations. The rotational distances and the 
translational distances in a given direction have used in [41] to formalize the constraints on the 
relative position of links in problems of modeling of mechanism's motion. 

The dynamic geometric situation is defined, for moving objects, by representing an objects 
as a four-dimensional sets in the space-time C; see [6], [39], and [41]. (Note that the distances 
between objects "in the space" and "in the time" are incomparable [64].) 

Let A* = UtG[o,i] [^(^)]' ^* = UtG[o,i] [Bit)] be the image of A, B in G\ We denote by ©* 
the partial addition "by the time", and by B* the reflection of B* with respect to the origin O in 
R^. Let next \K* be the cylinder in G"^ with a basis XK. Then the parametric family of objects 

ri{X,K,B*,A*) = U [A{t)®B{t)®XK] = A* ® B* (B XK* 
te[o,i] 

corresponds to the distance 71 A*) = inffg[o,i]{7i[i?(t), A{t)]} between A* and B*. 

It is clear that the suggested technique for solving the translational spatial planning problems 
can also be applied to geometric problems with the considered types of constaints on the relative 
position of objects. 

10 Computational issues 

From the observations of Sections 8 and 9 it follows that the C-space map of spatial plan- 
ning problem, is the region obtained by standard and/or regularized Boolean operations, and by 



26 



Minkowski operations on regular and/or non-regular objects. In general, the C-space map is a 
non-regular geometric object of various topology with non-manifold boundary; see Remark 5 
of Section 8. (Recall that the non-regular object may have external dangling faces/edges and/or 
isolated points, and/or internal entities such as cracks and/or isolated points [43].) Therefore 
for implementation of spatial planning we need the methods for representation and manipu- 
lation of such an objects. In this section we briefly consider the computer representations of 
geometric objects that are suitable for solving the spatial planning problems, and the strategies 
for computing the C-space maps. 

Representaions of geometric objects. The two representation schemes that are most widely 
used in solid modeling and computer graphics are boundary representation (BRep) and con- 
structive solid geometry (CSG) [43]. Let A be a point set of i?" (n = 2,3). CSG(A) is a 
Boolean composition of algebraic half spaces using regularized set operations. BRep(A) is a 
collection of closed faces/edges. The problems of CSG to BRep conversion and of BRep to 
CSG conversion have been studied in [43], [50], and [51]. 

The third type of representation scheme suitable for our purpose is the linear ray represen- 
tation (LRRep) [31], [39], denoted by LRR(A). (In [39] it is called the linear raster represen- 
tation.) LRR(A) is an approximation of an object A by a set of parallel segments belonging to 
a grid L of parallel lines, i.e., LRR(A)= An L. Conversions between BRep and LRRep, and 
between CSG and LRRep have been detally studied in [32]. 

The constructive non-regularized geometry (CNRG) methodology for representation and 
manipulation of non-homogeneous (i.e., made of several materials with different properties), 
non-closed point sets with internal structures and incomplete boundaries have been suggested 
in [48]. The work in [18] has proposed an approach for representation of geometric objects with 
non-manifold boundary. 

The boundary representation of non-regular geometric objects with non-manifold boundary 
using the techniques of [18] and [48] have been studied in [42]. It takes into account both the 
geometry and the topology of objects. The work in [42] have also considered the topological 
operations (complement, interior, closure, and regularization) on a single non-regular object. 

Boolean operations. Algorithms and implementation for computing the regularized set oper- 
ations on polyhedral objects have been proposed in [44]. Recall that the standard union AU B 
of two r-sets A and B always results in an r-set, whereas the standard intersection An B and 
the standard set difference A\B need not be regular: the set An B may have dangling edges, 
e.g., in case where A contacts with B along the portion of its boundary [43]; the set A\B may 
be partially open, e.g., in case where n i-B 7^ [36], [44]. Algorithms for computing the 
set operations on non-manifold boundary representation objects have been proposed in [18] 
and [46]. We assume below some familiarity with theory and algorithms of [18] and [44]. 

As mentioned in [44], the algorithms will work with curved objects, and they are insensitive 
to whether a solid's boundary is or is not a two-manifold, and is or is not connected. Hence, the 
algorithms of [44] can be modified to compute the standard Boolean operations on non-regular 
objects of various topology with non-manifold boundary. 

Let S = A0* B, where ® denotes one of the standard set operations. The main utilities used 
in algorithms of [44] to compute the boundary of S are the set membership classification (SMC) 
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[59], and the combining classifications, defined by means of the regularized set operations. 
See [44] and [59] for details. 

In [42] the algorithms of [44] have been modified for computing the standard Boolean 
set operations on non-regular (possibly unbounded) geometric objects of various topology, for 
BReps. To define and to combine the classifications the work in [42] have used the standard, 
but not a regularized set operations. 

Minkowski operations. Many various algorithms to compute the Minkowski operations have 
been proposed. Detailed surveys of previous work on computing the Minkowski operations can 
be found in [13], [20], [25], [26], [36], [60], and [63]. Algorithms for computing the Minkowski 
sums and the Minkowski differences in two and three dimensions are given, e.g., in [2], [3], 
[14], [15], [17], [25] - [29], [33], [37], and [60], for BReps, and in [31] and [39], for LRReps. 
Note that the referenced algorithms perform computing the Minkowski operations on regular 
objects. They can generate the manifold boundaries and are not applicable to the cases where 
the boundary of the resulting object is non-manifold; see [60] for details. 

In [12] and [20] have been presented an algorithms for robust and efficient construction of 
planar Minkowski sums for polygons using exact rational arithmetic. In contrast with most 
existing techniques the algorithms of [12], [20] directly handle the degenerate configurations, 
arising in the boundary of the Minkowski sum, such as internal isolated points and/or coinciding 
edges. In other words, these algorithms compute the outer envelope of A and B, i.e., the 
boundary of the open set iA(BiB (see subsection 2.1). The recent works in [61] and [62] have 
presented an algorithms for exact and efficient construction of Minkowski sums of polygons, 
and for exact and approximate construction of offset polygons, respectively, that handle the 
degenerate configurations also. Hence, the algorithms of [12], [20], [61], and [62] allow to 
construct the parametric families of polygonal objects used for computing the C-space maps. 

The algorithms of [12] and [20] are based on convex decomposition of polygons. However, 
as mentioned in [3], not all curved objects permit convex decomposition, e.g., an object with 
an inward concave edge. Therefore to handle the curved objects more suitable are the methods 
that deal with the geometric objects directly. 

In [42] the algorithms of [3] have been modified for computing the Minkowski operations 
on non-regular objects of various topology, for BReps, using the techniques of [18] and [48]. 

Distances between geometric objects. Many algorithms for computing the distances be- 
tween geometric objects have been developed (see [27]). Detailed surveys of previous work 
on computing the MTD between regular objects can be found in [10], [13], [16], [21], [27], 
and [35]. Algorithms for computing the distances concerning the outer relative position of 
objects in two and three dimensions (see subsection 2.2) are given, e.g., in [7], [10], [13] 
- [17], [35], [37], and [52], for BReps, and in [39] and [41], for LRReps. The algorithms 
of [12], [20], and [61] for robust construction of planar Minkowski sums can be used for com- 
puting the MTD between non-regular polygonal objects. 

Algorithms for computing the distances concerning the containment of objects in two and 
three dimensions (see subsection 5.1) are given in [39] and [41], for LRReps. Note that for 
this goal can also be used algorithms for computing the Minkowski difference, for BReps; see, 
e.g., [14] and [15]. Algorithms for computing the distances concerning the covering of objects 
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(see subsection 5.2) and algorithms for computing the distances concerning the containment of 
objects are similar. 

Algorithms for computing the translational distances in a given direction in two and three 
dimensions are given in [10], [13], [22], and [52], for BReps, and in [39] and [41], for LR- 
Reps, respectively. The algorithms of [39] and [41] are based on the partial vector operations 
(see subsection 9.1). Algorithms for computing the rotational distances and the partial vector 
operations in polar coordinates are given in [41]. 

Thus, using the algorithms for computing the various types of distances between geometric 
objects provide solving the generalized distance query problem, as defined in subsection 1.3. 
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